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Abstract 

We develop in the weak coupling approximation a quasi-non-Markovian master equation and study 
the phenomenon of decoherence during the operation of a controlled-not (CNOT) quantum gate in a 
quantum computer model formed by a linear chain of three nuclear spins system with second neighbor Ising 
interaction between them. We compare with the behavior of the Markovian counterpart for temperature 
different from zero (thermalization) and at zero temperature for low and high dissipation rates. At low 
dissipation there is a very small difference between Markovian and quasi no-Markovian at any temperature 
which is unlikely to be measured, and at high dissipation there is a difference which is likely to be measured 
at any temperature. 

PACS: 03.65.Yz, 03.65 Ud 

1 Introduction 

A quantum open system is generally characterized by a non unitary evolution of the reduced density matrix 
associated to of the central system and its interaction with the environment. Different types of approaches 
have been developed to understand the phenomenon of decoherence that arises in the open quantum systems 
which it is related to the lost of the interference terms of the product of the quantum wave function [T]- [5], 
that is, the non diagonal elements of the reduced density matrix go to zero value. Since the complete insulate 
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quantum system is almost impossible to have, decoherence becomes an intrinsic phenomenon related to the 
quantum principles and maybe related to the "emergent reality" of the classical world |14j-|18|. Many interest 
has been created in the phenomenon because of the difficulties it carries to perform quantum computation. 
Non-Markovian systems or systems where the environment is supposed to keep memory, is a topic in this 
subject and there is not a unified consensus about the best approach for studying the dynamics of this sys- 
tems [5]-|13|. which makes non-Markovian to be a very interesting subject. In the Markovian approach the 
non-unitary evolution equation is called "master equation" which is a differential equation for the traced over 
the environmental variables of the full density matrix. In principal, in the non-Markovian approach one will 
have to obtain an integro-differential equation for the density matrix and to establish the non-Markovian 
in it, but, how to measure non-Markovian ? it is still uncertain. Every approach needs to be intended to 
maintain the positiveness and trace equal to one for the reduced density matrix. The best known mathe- 
matical formalism which kept these conditions was given by Lindblad [1], who gave an abstract general non 
unitary evolution equation for the reduced density matrix, so keeping the Lindblad form in the equations is 
a good indication. We thought in quasi non-Markovian as an approximation to a master equation but with 
a temporal dependence in some of the coefhcients which defines the interaction with the environment which 
also depends on the Ising interaction between the spins and may lead to a different behavior of the traditional 
Markovian solutions. In addition, these solutions keep the completely positiveness of the density matrix. 

We use the weak coupling approximation for a system consisting of a linear chain of three paramagnetic 
atoms with nuclear spin one half [25^ , interacting with a thermal reservoir (not pure) consisting of a bosonic 
bath |19j-|21j. The temporal dependence in some terms, in the weak coupling approximation, is what we 
have considered as something beyond Markovian which is totally related to this type of system, and more 
specifically, to the Ising interaction between the nuclear spins. 

We study the decoherence of quantum controUed-not (CNOT) gates during operation in a quantum 
computer model. In this work, we are interested in determine the differences between the quasi-non-Markovian 
and Markovian behavior of a quantum controlled-not (CNOT) gate during its implementation on this model 
of quantum computer. In the first part of this work, we describe this model and the Hamiltonian of our 
quantum system interacting with a thermal reservoir, which consist of modes of an electromagnetic field in 
a cavity where the quantum system is. In the second part, we perform the weak coupling approximation to 
obtain a quasi-non-Markovian master equation. We want to point out that, even this model has not been 
built, it has been very useful for theoretical studies about implementation of quantum gates and quantum 
algorithms [251 HSl 130] which can be extrapolated to other solid state quantum computers [21] . Then, the 
analytical dynamical system of the reduced density matrix elements are obtained, and the results of the 
numerical simulations are presented. We present mainly the differences between the Markovian and the 
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quasi-non-Markovian behavior on the reduced density matrix elements. 



2 Hamiltonian of the system 

The Hamiltonian that describes the ideal insulated system of a linear chain of N paramagnetic atoms with 
nuclear spin one half inside a magnetic field 

B{z, t) = (6cos(wt + (p),bsm{ujt + ifi),B{z)) (1) 

where b, uj, and (p are the amplitude, the angular frequency and the phase of the RF-field, and B{z) represents 
the amplitude of the z-componcnt of the magnetic field, is given by [2S] 

TV N-1 N-2 

Hs ^ - ^flk ■'Bk + J ^ S^S^+i + J' SlSl^2 > (2) 

k=l k=l k=l 

where flk represent the magnetic moment of the kth-nucleus, which it is given in terms of the nuclear spin as 
flk = 7(5'fc J S^), with 7 being the proton gyromagnetic ratio and Sf. being the jth-componcnt of the spin 
operator, B^, represents the magnetic field Eq. ([l]) valuated at the location of the kth-nuclear spin (z = Zk). 
The parameters J and J' represent the coupling constant at first and second neighbor interaction. The angle 
between the linear chain and the z-component of the magnetic field is chosen as cos 9 = l/\/3 to eliminate 
the dipole-dipole interaction between the spins. 

We can write the Hamiltonian ^ in its diagonal and non diagonal with respect a chosen basis in the 
z~ projection as 

Hs = Hq + Hrf (3) 



where 



and 



N N-1 N-2 



k=l k=l k=l 



n 



N 



k=l 

Here we have that: ujk = ^B(zk) is the Larmor frequency of the kth-spin, Q — jb is the Rabi frequency, 
and — ± represents the ascend operator (+) or the descend operator (— ). The Hamiltonian Hq is 
diagonal in the basis {\aN ■ ■ ■ c^i)} with = 0, 1 (one for the ground state and zero for the exited state ). The 
action of the spin operators on its respective qubit is given by S^\ak) = l)"*" |afe)/2, S'^\ak) — hda^.^i\0), 
and S'i^\ak) — hSa^flll)- The eigenvalues of Hg in this basis are given by 

f- ( N N-l N-2 ~j 

= -2 E(-i)"''^fc + E (-1)"''+"^+^ + J' E (-1)"'^+"'"+^ . (6) 

U-=i fe=i fe=i J 
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The elements of this basis forms a register of N-qubits with a total number of 2 registers, which is the 
dimensionality of our Hilbert space. The allowed transition of one state to another one is gotten by choosing 
the angular frequency of the RF-field, w, as the associated angular frequency due to the energy difference of 
these two levels, and by choosing the normalized evolution time fit with the proper time duration (so called 
RF-field pulse). The set of selected pulses defines the quantum gates or the quantum algorithms, and CNOT 
quantum gate is the gate we want to study. 



Consider now this system to be immerse in a " mixed thermal bath of oscillators " such that the Hamiltonian 
of the bath is of the form 

oo 

He = ^/iWja]aj. (7) 

j 

The Hamiltonian of the interaction between the central system and the bath will be taken in the form 

Hi = Y^ au,S,E, = ^ (a\=Sta, + a\=S-a]) , (8) 

kj k,i 

where the operator Ej is defined as Ej = aj +a], Sk is the polarization operator, Sk = + , and we have 
taken into account the Jaymes-Cummings rotating wave approximation for the interaction [22] . in order to 
considerer an excitation-de excitation process of the system trough the coupling with the bath of oscillators 
with characteristic frequencies near the resonant frequencies of the transitions. The constants a^'' , i — 1...3 
are phenomenological parameters that measures the coupling between the system and the environment and 
aj(a|) are the rising (lowering) operators in jth number of photons in the bath. We can write the total 
Hamiltonian as 

H = Hd+Wi , (9) 

where and Wj are given by 

Hd=Ho + He (10a) 



and 

Wj{t) = Hrf + Hj = --^Y. {^^^'^''"'^3+ + e-("*+'^)^-) + ^ {a1^S+a, + 4^S-a]) . (10b) 

fe = l k,j 

3 The weak coupling approximation 

Now, for dealing with the non ideal situation we start with the dynamical equation of the evolution of the 
density matrix for an initially decoupled state in the system plus the environment 

p^ps®(7E (11) 
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where ps is a pure state of the central system and ge is a thermal stationary mixed state of the environment. 
In the interaction picture the equation of evolution for the reduced system is 

j/s{t) = -'-TrE{[Wi{t)rps{t)®aE\}, (12) 
where in this interaction picture one has 

ps[t) = e'"°''^pse-'"°''^ , (13) 

and 

~a,[t) = e'^^^l^a^er'^^^l^, ~a\{t) = e^^^l^a\e-''^^''^^ (14b) 
with the operator 51^ being defined as 

=u^- J{Sl^, + S^) - J'(5,%2 + Sl_,), (15) 

which commutes with the Hamiltonian Hq. The eigenvalues of this operator Qk, 

nk\t) , (16) 

are given in the appendix. 

The time integration of the system in the interval [t, t + At] is given as follows 

■ rt+At 

psit + At)^ps{t)-- dhTrE{[Wi{h),p{h)<E,<jE]}, (17) 



Then by doing a successive change of variables and substituting in (17), up to second order terms, using 



Markov approximation and Eq. (11), we obtain 



2 rt+At / 1 \ ^ rt+At rti 

Aps = ^ dhTrE{[Wi{h),ps{t)®aE]}+i^j / dh dt2TrE[Wiih),[Wiit2), psi^) ® '^e]], 

(18) 

where time locality is shown inside the integration with the term Ps(t)j s^nd we have set Apg ~ ps{t + At) — 
Ps{t)- One would expect that within this weak coupling approximation, the interaction of the central system 
with the environment would show a perturbation to the closed system. By substituting the corresponding 



time dependence form of Wi{t) in ( 18 ), one can sees that the following relation must be satisfied (notice that 
rijj*^ « ujk for all basic state \i)) 

\iu + n^^^\At «l, fori=l,...,8 (19) 

which determine the time path length where there is no interaction with a time dependent external field and 
no interaction between the spins How smaller this path has to be is not resolved, but definitively not that 
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small compared to the relaxation times of the environment te such that the Markov approximation still being 
valid. The lost of the separability of the initial system-environment states Ps® pE for a smaller Ai could exist 
since a longer time will have to pass for the evolution in the system and therefore correlations between the 
system and environment can arise, but in the case when we have a thermal state for the environment which 
is our case, any correlation generated by the evolution of the central system will rapidly decay. Integrating 



(18) and imdcr the condition (191, the master equation takes the form 



^PS^^) - ^._[H,f{t)rps{t)] + ^,(^ I dT r dhTrE[Hiih),[Hi{h-T),psit)(g>aE]], (20) 



At ih' " At\ih 



Jt 



where we have made the change of variables t2 = ti — t with t £ (0, At) such that ti € {t + T,t + At) and 



divided all by At. The first term in the right hand side of (20 1 describes the ideal part of the dynamics in 



the interaction picture (von Neuman evolution), and the second part describes the open dynamics. 

For a thermalized mixed environmental system one can sees that {a(s)a{t))E = {a' {s)a^{t))E = 0, then 
by doing typical calculations consisting in integrating over ti by using the spectral representation of the flk, 
performing the wave rotating approximation and regrouping terms, it follows that 

^ = l^[Hrf{t),ps{t)] -^Y.fdT {(A,(r)4)^e-"'^- {S^S^psit) - 5,7ps(t)^,t) 

+ (4(T)Afe)Be^"'=- {S-S+ps{t) - S+ps{t)S-) + h.c.} (21) 

where the limit At — > dt has been taken, the superior limit in the integrals has been put infinity since the 
correlation functions decay exponentially in time, and the following definitions have been made 

Ak^J2 9^ka^, 4 = 9*ik4 ■ (22) 

The coefficients gik and g*^. are related to the coupling of the central system with the environment and 
depends on the characteristic frequencies of the modes in the neighborhood of each spin. The correlation 
functions are described by the Fourier transform of certain spectral density, j{uj), associated to the continuous 
modes in the thermal bath, 

/■OO 

{A{t)A'^)e = 7o / dLoj{u)e'"^. (23) 



(24) 



with 7o = Igp. The correlation functions appearing in the a factor can be written as 

dTe^'''''^A,{±T)Al)E = -Ik ± ifk. dTe±*"'=-(4(±r)Afe)i. = -ll T *f ^ , 

where we get the operators 
7fc = 2^7„j(^^fe), 7i = 2^7o/(f^fc), ffc=7„P.V. / T^r^d^, f ^ = 70P.V. / /^d^, (25) 



6 



being P.V. the Cauchy principal value. These operators are diagonal on the above basis, 7fc|j) = 7fc'^l*) for 
example, and their eigenvalues are denoted with a upper index (see appendix). 



By regrouping terms in (21 ) and going back to Schrodingers picture, we obtain the following master equation 



j^Ps^l^[Hs + HL,ps]+Cp{t), (26) 



with Cp{t)i defined as 



7fe 



fe 



t 



2 



;,t 



PS 



S^Si -T.^t''\t)St\n)p's-'-\MS^ ^ (27) 



where p^"'™-* = {n\ps\m) are the matrix elements of the initial reduced density matrix operator, and in 



Eq. ( 26 ) is given by 

HL^hS^S+-riS+S^, (28) 

which represents a Lamb shift Hamiltonian and can be not considered in the dynamics since it commutes 
with the entire Hq of the central system. In addition, it only generates a global shift in the spectrum. The 
time dependent coefficients are explicitly given by 

7fc (t) = e V ^ ; , (29) 

and they represent local phases for the non diagonal terms of the equations of the density matrix. Therefore, 
the positiveness and trace equal 1 are still satisfied for the density matrix. These phases depend linearly on 
the Ising coupling constants and bring about the quasi non-Markovian behavior of the system. 

If we consider low Ising coupling with respect the Larmor frequencies, then we can make the following 
approximation 

(30) 



and (26) takes the following form 

^^Ps = ^^[Hs.Ps]+Cp. (31) 
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where the term Cp2 is given by 



^-"2 = 7^ E { y {S+S^ps - 2SlpsSt + PsStS^) + 1^ {S^S+ps - 'iS+PsS^ + ps^.75+) | 



(32) 



with the coefficients 7^ and written as 



7fc = 2TTjoj{uJk), ll = 27r7oj^(wfe). (33) 

This type of master equations generates no-correlated thermaUzed cases which describes spontaneous and 
thermaUy induced emission-absorption process [3], [13], [231 ■ The environment generates excitations or de- 
excitations in the closed system by absorbing-emitting photons of the thermal bath. In this work, we want 



to establish the differences between Eq. ( 26 1 which may describe a quasi- non-Markovian process via the 



oscillating term in the non diagonal elements of the dissipator, and Eq. (31 1, which is the typical Markovian 
master equation for a system immerse in a bosonic field. 



4 Physical quantities 

Let us considered a thermal bath of radiation modes at a temperature T. The environmental density matrix 
is given by 



The interaction Hamiltonian between the central system and the environment is represented by a coupling 
between the polarization operator and a bosonic modes operators. The correlation functions involved in the 
system ((ylfc(±T)ylj.)£;, {A!^^.{±t) Ak) e) are calculated, 

{A,{±r)A\)E = E \9\%e^'-^^ ( .^^Ii;;^ + A ^ (4(±-)^.)i. = E l^lf. ^^^Z^^TT^I- (^5) 

The sum over i is dense (there are an uncountable number of radiation modes). If the volume containing this 
modes is large enough, we can go from a discrete distribution to a continuous distribution of the characteristic 
frequencies of the radiation modes. The number of characteristic frequencies with wave vector components 
/ in the interval dfxdfydfz in the volume V is given by 

. ^d., (36) 

where / = c-lu. Thus the sum in the correlation functions can be changed by an integration over de frequencies 
with the proper weight factor, 

{Ak{±T)Al)E - ^ r du;c.3(7V(^) + i)eT--, {AI{±t)A,)e = ^ H du:uj^ N {u)e^'^^- (37) 
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where we have taken a hnearly dependence on the characteristic frequencies of the radiation modes, \g\ 
If/fcpWj, and the Planck's distribution function, 



jk 



N{uj) 



1 



(38) 



Comparing this resuhs with the definitions in ( 24 ) and ( 25 ) , one can sees that 



(39) 



and 



7o 



V\g\l 



(40) 



Once we get the definitions of all these constants, we can proceed to solve the above equations. The evolution 
equations of the matrix elements are given in appendix. 



5 Simulations and results 

Our registers are made up of three qubits \ ABC) with A, B,C = 0, 1, or written them with decimal notation, 
|1) = |000), |2) — |001) and so on. The parameters used for our simulation are taken from i30j regarding 
the Larmor frequencies of the nuclear spins and we take a higher Ising coupling strength for modelling the 
differences between the Markovian and the quasi-non Markovian regime but maintaining the 2tt k method 
[3D]. These parameters are (in units of 2tt MHz) as 

ujA = 400 , OJB = 200 , = 100 , J = 25 , and J' = 1 (41) 

There is still one free parameter which is the strength of the coupling between the environment and the central 
system This will allow us to model high or low dissipation rates of an homogeneous or inhomogeneous 
environments. We take the assumption that the environment is acting homogeneously on each qubit, that is, 
there is a set of baths of characteristic frequencies affecting more closely the resonant frequencies of each spin. 



The reduced density matrix is then made up of 8 x 8 complex elements, and if the initial state is always taken 
as the exited state |1) — |000), this means that the initial reduced density matrix has the values pn — 1 and 
Pij = for i,j ^ 1. 

5.0.1 Controlled-Not (CNOT) quantum gate 

To get the CNOT quantum gate starting from the ground state |1) = |000), one applies a tt/2- pulse between 
this state and the state |3) = |010), with resonant frequency uj — ujb — J, to get the superposition state 
(|1) + |3))/-\/2. Then, one applies a resonant 7r-pulse between the states |3) and |4) — |011), with resonant 
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frequency uj = ojc + J/'2, — J'/2, to get the final desired state (|1) + |4))/\/2 which means that the expected 
CNOT density matrix would be such that pn = pi/^ = p4i = p^i = 1/2, and all the other elements are equal 
to zero. In addition, one allows the system to have two and a half more resonant 7r-pulses to have a better 
look of the CNOT behavior. 



5.0.2 Dynamics at room temperature 

We start modeling the dynamics by considering that the environment is at room temperature (T = 300 
Kelvins) . This assumption will make the system to evolve into a thermalized mixed state. We present in the 
following figures the differences of the behavior of the diagonal terms and the coherent terms involved in the 
CNOT quantum gate. 

Figure 1 shows the behavior of the diagonal elements of the CNOT gate for low (7 ft! J' x lO"*^) and 
high (7 « J' X 10"-'^) dissipation rates (the high dissipation rate is still in the limit of considering the 
approximation of a perturbation of the central system) for the Markovian and semi(quasi)-non Markovian 
regimes at a temperature T = 3Q0K. We can see that the differences between them is very small but still 
distinguish. For each case, Markovian and quasi-non Markovian, the environment will lead the central system 
into a thermalized mixture states, with a thermalization time depending on the coupling constant with the 
environment. We need to point out that this difference increases as the spin coupling constant a first neighbor 
increases its value. 
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Diagonal elements of density matrix for CNOT quantum gate 




Diagonal elements of matrix density for CNOT quantum gate 

Markovian regime for high dissipation at T-300 Kelvins 




Diagonal elements of matrix density for CNOT quantum gate 

Semi-non Markovian regime for low dissipation at T— 300 Kelvins 




Diagonal elements of matrix density for CNOT quantum gate 

Semi-non Markovian regime for high dissipation at T-300 Kelvins 

-r 




FiGURE 1 Diagonal elements of the density matrix for the CNOT quantum gate for low (left) and high 
(right) rates of dissipation in the Markovian and quasi-non Markovian regime at T=300 Kelvins. 



Figure 2 shows the behavior of the coherent terms involved in the CNOT gate. We can see that for high 
dissipation rates a fast thermalization of the system, making the coherent terms goes to zero very rapidly. 
The term (Ipjia) is related to the first pulse which makes the superposition state for the CNOT formation. 
Therefore, it has higher amplitude, since it will take some time for the environment to completely thermalize 
the whole system. For the last 7r-pulse for the CNOT formation, we see that the decoherence is already high. 
There is a very small sudden birth of coherence in the term |p|i4, due to the pulses of the magnetic field 
needed to perform the quantum gate. For the low dissipation cases (7 « J' x 10^^), the quasi-non Markovian 
and the Markovian regime are very similar, yet we can see a very small differences in the amplitude of the 
coherent terms. The elements involving the higher energy level ( |1) = 1 000)) their amplitude seem to have 
a grater amplitude for the quasi-non Markovian regime than in the Markovian regime. This situation is 
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contrary in the element |p34|. 



Coherent terms lp|^^ 

low dissipation al T=300 Kclvir 



Semi -non Maikoviau 





Coherent terms Iplj^ 

high dissipation al T=300 Kelvins 
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low dissipation at T=300 Kelvins 




Coherent terms lp|^^ 

high dissipation ai T=300 Kelvins 
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Coherent terms Ipl^^ 

low dissipaUon al T=300 Kelvi 




Coherent terms lp|^^ 

high dissipation for T=300 Kelvins 
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Figure 2 Coherent elements Ipial, |pi4| and \p34\ of the density matrix for the CNOT quantum gate for 
low (left) and high (right) rates of dissipation in the Markovian and quasi-non Markovian regime at T— 300 
Kelvins. 
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5.0.3 Dynamics at the thermal vacuum 

At T = Kelvins, the master equation takes the form as described in |23| for the A-Independent environment 
cases, and we still have the time dependent terms on the non diagonal elements of the master equation, 
referring to the quasi-non Markovian case. Figure 3 shows the behavior of the diagonal elements of the 
CNOT quantum gate at T = Kelvins. We can not see a distinguished difference between the Markovian 
and quasi-non Markovian regimes as we did in Figure 1. However, at high dissipation rate, we still seen for 
both cases (Markovian and quasi-non Markovian), the rise of the equilibrium ground state at the end of the 
whole process. This happens because our initial state is the most exited state, and during the process of 
dissipation, the environment is not giving off any energy to the system, the quantum system will deliver the 
energy to all the other states, exiting them. Therefore, by dissipation, all of them go back to zero, leaving 
the system in the ground state |8) = (purple curve in the figures). 



Diagonal elements of matrix density for CNOT quantum gate 
Diagonal elements of matrix density for CNOT quantum gate Markovian regime tor high dusipauon ai t=o KeMns 




Figure 3 Diagonal elements of the density matrix for the CNOT quantum gate for low (left) and high 
(right) rates of dissipation in the Markovian and quasi-non Markovian regime at T=0 Kelvins. 
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Figure 4 shows the behavior of the coherent terms of the reduced density matrix at T = O^iT. For 
low dissipation rates (7 ~ J' x 10^^), we see a similar behavior as in the room temperature cases. For 
high dissipation rates ^ J' x lO^-"^), we see more significant differences between the Markovian and the 
quasi-non Markovian cases, suggesting that this effect could be observable experimentally. 



Coherent terms lp|^^ 

low dissipalioii al T=0 Kelvins 




Coherent terms lp|j^ 

high dissipalion al T=0 Kelvins 



10 20 30 40 50 60 70 




Coherent terms lp|^^ 

low dissipation at T=0 Kelvins 



Coherent terms Ipl^^ 

high dissipation at T=0 Kelvins 





Coherent terms lp|^^ 

low dissipation at T=0 Kelvins 



Coherent terms lp|^^ 

high dissipation at r=0 Kelvin 





Figure 4 Coherent elements Ipial, \pia\ and |p34| of the reduced density matrix for low (left) and high 
(right) rates of dissipation in the Markovian and quasi-non Markovian regime at T = 0°K. 
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5.0.4 Purity calculations. 

The purity function, P{t) = tr(p'^), is a measure of Irow close a quantum system is from its description as 
a pure state quantum system (the density matrix be written in term of a wave function p = and 
varies between 1 and 1/d (d the dimensionality of the density matrix) . This function may decay with the 
decoherence since the system may move away from an initial pure state. Therefore, this function can be used 
to characterize the environment. 

Figure 5 shows the behavior of the purity for the CNOT gate at room temperature and at T = 0°K. At 
room temperature it is observed a thermalization of the system, and at T = 0°K a recovery of the purity 
since the system goes to the quantum ground state, depending on the dissipation rate. 

Purity of CNOT quantum gate Purity of CNOT quantum gate 

low dissipation at T-300 Kelvins high dissipation at T-300 Kelvins 




i=xia l=x/£J 



Purity of CNOT gate Purity of CNOT quantum gate 



low dissipation at T=0 Kelvins high dissipation at T=0 Kelvins 




t=T/Q l=z/Q 



Figure 5 Purity for the Markovian and quasi-non Markovian regimes at T=300 and T=0 Kelvins for low 
(left) and hight (right) dissipation rates. 
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6 Conclusions 



Within the weak coupUng approximation for the study of quantum discrete system with environment, we 
have obtained a quantum master equation with a time dependent non diagonal dissipative coefficients which 
shows a quasi-non Markovian behavior. We have solved numerically the master equation for the reduced 
density matrix associated to our linear chain of three nuclear spin system interacting with the environment. 
We have made the simulation of CNOT quantum gate operating in this dissipative environment and within 
the validity of this approximation. We have study the behavior of system-environment interaction with 
this quasi-non Markovian master equation and compared the results with the Markovian counterpart. The 
decoherence of this quantum logic gate have been determined, and we have seen a different behavior of the 
decoherence with the quasi-non Markovian and with Markovian master equations. 

This difference between quasi-non Markovian and Markovian approaches on the reduced density matrix 
elements grows with the dissipation coefficients defined in the master equations, but the diagonal elements 
remain almost identical over the two types of process. Therefore, for low dissipation the measuring apparatus 
will not bring any information of the environmental interaction for a Markovian or quasi-non Markovian 
process, and for high dissipation this difference can, in principle, be measure. In addition, this difference 
must increase as the spin coupling parameter a fist neighbor increases since the spectrum becomes much 
more well defined. This comparison was also made using the purity parameter. For strong dissipation at 
T = Kelvins, we found that purity may increase because, the condition trp = 1 on the density matrix, and 
this implies excitation of the equilibrium state involved in the dynamics (ground state), causing the system 
to try to return to a pure quantum state description. 
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7 Appendix 



The evolution equation for the density matrix elements are given from Eq. ( 26 1 by 

'^-P^^-Uh,pU+\c{p)] , a,/3 = l,...,8. 
at n r. 

Making the following definition 

(v^)q/j = ^[-ff,p]Q/3 and Cp^^ = ^ £{p) 

one gets 



a/3 



(Xl) 



{X2) 



Von Neuman (vN) part. 



(vN),, = -^e*M+'^) (p,, + + + (pi2 + Pi3 + Pis) 



(Vl) 



(vN) 



= - (^c - j/2 - j/2) - ^e^(-*+^) (p22 + P32 + P52 ' Pii) + ^e-»("*+'^) (pi4 + Pie) (^2) 
(i;7V)i3 = - (^B - i) Pi3 - (P23 + P33 + P53 " Pii) + ^e-'("*+^) (pi4 + pit) (^3) 



{vN),^ ^^(^ujB+oJc- j/2 - //2) pi4 - ^e^(-*+^) (p24 + P34 + P54 " Pi2 ~ P13) + ^e-'("*+^Vi8 (V^4) 
(^^A^)i5 = - (c^A - i/2 - j'/2) - ^e^('^*+'^) (P25 + P35 + P55 - Pii) + ^e-^("*+^) (pi7 + Pie) (^5) 
{vNU = -{u;A+iOc- j) Pie - ^e'^^'+^'> (pse + P3e + Pse - P12 - P15) + ^e-'("*+^Vi8 {V6) 



n 



{vN),, = -{u:a+ojb- j/2 - //2) p^j - ^e^("*+^) {p,, + p,, + p,, - p,, - p,,) + ^e-*('-*+'^)pi8 {V7) 

m 

{V9) 



{vN)is = - {UJA + UJB + ^c) Pis - —6*^'^*+'^) {p2S + P38 + P58 " Pl4 ^ Pl6 ^ Pi?) 



{vN)22 = -^e'("*+'^) (p42 + Pe2 - P2i) + ^e-*('^*+'^) {-pi2 + P2e + P24) 



{VN)23 ^-{WB-UJC- j/2 + 3/2^ P23 - ^e*(-*+'^) (p43 + Pe3 - P21) + ^e-*('^*+'^) (-P13 + P24 + P27) 

{VIO) 

{VN)24 = -0JBP24 - |e^("*+'^) (p44 + Pe4 - P22 - P23) + |e-'('^*+^) (-P14 + P28) (Vll) 
{vN)25 = -{iOA- UJC) P25 - ^e^('"*+'^) (P45 + P65 " P21) + ^e-'("*+'^) (-pig + P27 + P2e) (1^12) 



{vN) 
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(loa - j/2 + 3 12) P2e - ^e^("*+^) (p46 + Pee - P22 - P25) + ^e'^^^'^^^ (-pig + P28) (^^13) 

(VN)2, = ~{UJA^^B- CJc) P27 ~ |e^("*+^) (p47 + Pe7 - P23 ^ P25) + ^e-'("*+^) + P28) (^^14) 



(uJA + C^B + j72 + 3 12) P28 - ^e'(-*+'^) (p48 + P68 - P24 - P2e - P27) + ^e-^("*+^) (-P18) 
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(yi5) 
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(^^A^)33 = --e*('^*+'^) (P43 + P73 - P3l) + -6-'^^'+^^ (-pi3 + P34 + P37) (1^16) 
{VNU = - {oJc+j/2 - 3 12) P34 - |e^(-*+'^) (p44 + P74 - P32 " P33) + ^^-^^'^^^ (-pi4 + P38) (V'lT) 

(^^A^)35 = - - a;B + i/2 - j'/s) P35 - ^e^(-*+^) (^45 + P75 - P31) + ^e-*(-*+^) (-pis + P37 + P36) 

(F18) 

(wAr)36 = - (W^ - + Wc) P36 - ^e*("*+'^) (P46 + P76 - P32 - P35) + ^e-*("*+^) (-pi6 + P38) (1^19) 

{vN)i^ = - + j/2 - j72) P37 - |e^("*+^) (^47 + P77 - P33 - P35) + ^e-*(-*+'^) (-pir + Pas) 

(^;iV)38 = - (WA + + j) P38 - ^e'('^*+'^) (P48 + P78 - P34 - P36 - P37) + ^e-^('^*+^) (-pig) (1^21) 

(^;iV)44 = -^e^('^'+'^) (p84 - P42 - P43) + ^e-*("*+'^) (-P24 " P34 + P48) (^23) 
(«Ar)45 = -{UA-'^B- ^C) P45 - ^e'("*+^) (P85 - P4l) + ^e-'(-*+'^) (-p25 - P35 + P47 + P46) (1^24) 
(t;iV)46 = - - - j/2 + j72) P46 - ^e^('^*+^) (pse - P42 - P45) + ^e-*(-*+'^) (-p26 - P36 + P48) 

(y25) 

(^;7V)47 = - (u;^ - wc) P47 " ^e'("*+^) (p87 " P43 " P45) + ^e-^("*+'^) (-P27 " P37 + P48) (^^26) 
{VN)^^ = - (WA + i/2 + i72) P48 - ^e^(-*+^) (P88 - P44 - P46 - P47) + ^e-^(-*+'^) (-P28 - P38) (V^27) 

(^iV)55 = -^e'(-*+'^) (p65 + P75 - P5i) + ^e-*(-*+'^) (-pis + P57 + Pse) (^28) 

= - (WC - j72 + i72) P56 - ^e^(-*+'^) (P66 + P76 - P52 - PSb) + ^e-^(-*+^) (-P16 + P58) (1^29) 

(t;7V)57 = -a;BP57 - ^e*(-*+'^) (p67 + P77 - P53 " Pss) + ^e-^('^*+'^) (-pir + Pss) (^30) 

(t;Ar)58 = - (a;s + + i/2 + i72) P58 - ^e^("*+'^) (p68 + P78 - P54 - P56 - P57) + ^e-'(-*+'^) (-P18) 

(1/31) 

(^^A^)66 = -^e^(-*+'^) (p86 - P62 - P65) + |e-^(-*+'^) (-p26 - P56 + Pes) (1^32) 

(t;iV)67 = - (wB - '^c + i/2 - 3 li) P67 - ^e*('^*+'^) (p87 - P63 - Pes) + |e-^("*+'^) (-p27 - P57 + Pes) 

(^33) 

{vN\^ = - (WB + j) Pe8 - ^e^("*+^) (p88 - Pe4 - Pee - Pe7) + ^e-^('-*+'^) (-p28 - Pss) (^34) 

{VN),^ = -^e^(-*+'^) (P87 - P73 - P75) + ^e-^(-*+'^) (-P37 - P57 + P78) (1^35) 
(t;iV)78 = - (wc + j72 +i72) P78 - ^e^(-*+'^) (p88 - P74 - P7e - P77) + ^e-^("*+^) (-P38 - P58) (1^36) 

(^;iV)88 = -^e'("*+^) (-P84 - Pse - P87) + |e-'('*'*+*') (-P48 - Pes - P78) {V37) 
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Dissipation part. 

i-Pl2 - - I ^ \ ^ \ ^ I Pl2 H ^ Ja \t)P56 + ^ 7b (*)P34 

/^/l) _L (1) t(3) ^(1) , ^,(3) \ t(l) , t(3) t(l) , t(3) 

-t-PM - - I ^ \ ^ h ^ I Pl4 H ^ 7a Wp58 

/Jl)_L^t(5) ^(1) ,^(5) ^(1), ^(5)\ ^t(l) , -,t(5) ^ ^ ^t(l)_L^t(6) ^ ^ 

i-pie - - I ^ ^ ^ 1 ^ I P16 H ^ 7b (*;P38 

i-Pi? = - I ^ ^ ^ H ^ I Pl7 H ^ 7c (*)P28 

rpi8 - - ^ + ^ + ^ J P18 

i-P23 - - I ^ \ ^ 1 ^ 1 P23 H ^ 7a (^jP67 

,^(4) ^(2) , -,t(4) -,t(2) , t(4)\ ^t(2) , t(4) ^t(2) , ^(4) ^ ^ 

- I 7a +7a i 7b +7b 7c +7c \ „ , 7a +7a „, 86 r^^.^ , 7c +7c i3 r^'.^ 

i-P24 - - I ^ 1 ^ ^ ^ I P24 H ^ 7a l*JP68 H ^ 7c Kt)Pl3 

/-n - I 7a +7a I /B + IB I /C + \ , I /B + 7_B „ 74W V 
''-P25 - - I ^ ^ ^ ^ ^ I P25 H ^ 7b irjP47 

/^(2) t(6) ^(2) , (6) ^t(2) , t(6)\ ^t(2) ,^t(6) , ^ ^(2) , ^(6) , , 

rp27 - - ^ + ^ + ^ J P27 

/ (2) , t(8) (2) , t(8) t(2) , t(8)\ (2) , (8) 

rn - ( 7a^ +7a I 7b' +7b , 7c +7c ] ^ , 7c + 7c .,(17)., . 

''-P28 - - I ^ \ ^ H ^ I P28 H ^ 7a WPi7 

r ( (3) , t(3) , (3)\ , t(3) , (3) , t(3) 

-Cpss ^ -yiA +7s +7c j P33 +7a V77 + 7b Pll +7c P44 

(4) ^,t(3) , t(4) ^(3) , t(4)\ ^t(3) , t(4) ^,t(3) , t(4) ^ ^ 

- I 7a + 7a , 7b + 7b , 7c + 7c \ ^ , 7a +7a ,, 87 ^^ j. 7s_jt7^^, 12 c .x 
^-^34 = - ^ H ^ \ P34H ^ 7a WP78H 7. 7s WP12 
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/ (3) , t(5) t(3) , (5) , (5)\ t(3) , t(5) 

-^Pss = - I ^ \ ^ 1 ^ I P35 H ^ 74 ( 



+ 7c (64) . , N 
2 7a (*)P46 



/ (3) , t(6) t(3) , (6) (3) , t(6)\ 

r„ - { 7a +7a i 7b ' ^Tb i 7c' +7c \ 
'^P'' - - \ 2 + 2 + 2 J 



P36 



2 2 2 



\ ^(3) (7) ^t(3) t(r) 

1 P37 + ^^^7i^^Hi)Pl5 + ^^^^ 



7c^H*)P48 



'38 = — 



^(3),^t(8) ^t(3) t(8) ^(3) , t(8)\ ^t(3) t(8) 

^ 7b +7b ^ 7c +7c ) ^ 1b +1b ^,(i6)f 



P3S+ 7rft)pi6 



/ (4) , t(4) , t(4)\ , t(4) , (4) , (4) 

Cpii = - \^^\> + '^'^ > + ^ j P44 + 1a P8S + Tb P22 + 7c V33 



/^W + ^t(5) t(4) (5) t(4) (5) ■ 



P45 



'46 



|^ 7in^7r , 7r+7i^^ , 7^+7^^^^ ' 



P46 



(4) , (6) 
7c +7c ./35) 



7r(i)P: 



'35 



(4) , (7) 
n I 7b +7b .X25)/,n^ 
P47 H ^ 7b WP25 



%(4)+^t(8) ^t(4) , t(8) ^t(4) t(8)\ ^(4) 1^(8) ,,,, -Y^^^ + -y^^^^ 

7a +7a I 7b 7b I 7c + 7c | „ , 7b +7b „(26^^\„ i 7c_Jl7c^ 37 . x 
^ \ ^ 1 ^ /348 H ^ 7b WP26 H ^ 7c WP 



'37 



rp55 = - (7^ + 7f ^ + 7^^>55 + 7iV + 



t(5) , t(5) 

7b V77 + 7c pee 



7a +7a I 7b + 7b i 7c +7c \ . 7a + 7a ^12). -. . 7b + 7b ^, 87 / . 
^ 1- 2 2 j 2 ^ -^^^^ 2 ^^^'^^ 



+ 7c' \ ^ , 7a ' + 7a ^ ^.(13) . ,^ ^ , 7c" +7c'\ .(86) / .^ ^ 

2 J ^^'^ 2 ^^^^^ 2 ^^^''^ 



2 2 2 

/^t(5) , t(8) ^(5) , t(8) ^(5) , t(8)\ -^(5) (8) 

— / 7a +7a I 7b + 7b I 7c + 7c \ ^ , 7a +7a ^, i4 /.y^ 
i-P58 - - I ^ 1 2 ^ 2 I 2 

r ( t(6) , (6) , t(6)\ , (6) , t(6) , (6) 

-Cpee = - (^7a + 7b + 7c J Pee + 7a V22 + 7b V88 + 7c Vss 



'67 



/^t(6) , t(7) (6) t(7) ^(6) , t(7)\ ^(6) +^(7) 

: - ' + V"" + j P67 + 7f^(*)P: 



2 ' 2 ' 2 

t(6) t(8) ^t(6) , t(8) t(6) , t(8) \ (6) (8) (6) (8) 

rp77 = - (7f ^ + 7^^^^ + 7^0 P77 + 7i^V33 + 7i^55 + ' 

7r+7r , 7^+7^ , 7^^^+7r ' 

2 2 " 



t(7) 
■7c V88 



^ 1 P78 H 2 V-)P34 H 2 V-)P56 

\ , (8) , (8) , (8) 

) P88 + 7a V44 + 7b V66 + 7c V77 



r ( t(8) , t(8) , t(8)\ , (8) , (8) , (8) 

'Cp88 = - \1a +7b +7c ) "-^-'^ -^-^ ' 
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Eigenvalues of the Vlk operator. 

The eigenvalues equation is written as = for a three nuclear spins \ABC). The basis is taken 

in decimal notation, like |1) = |000), |2) = |001), and so on. 

l^W ^ ^(5) ^ _ J/2 - j'/2, = 17^^ =ojB-J, ^''c = ^c' = wc - J/2 - f/2, {Al) 

=flf =LJA-J/2 + J'/2, 0(3^) =0g) =a;B, = q'^^^ = coc + J/2 - J' /2, {A2) 

= n^l^ =ujA + J/2- J' /2, n^g'' ^n^p =ujB: n^-^^^ = np =ujc - J/2 + j'/2, (as) 

fiW ^ ^(s) =u;a + J/2 + J'/2, = nf =UJB + J, ^'c = = + J/2 + J'/2, (A4) 
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